Removing motion and physiological artifacts from intrinsic BOLD fluctuations using short echo data

AbstractDiffering noise variance across study populations has been shown to cause artifactual group differences in functional connectivity measures. In this study, we investigate the use of short echo time functional MRI data to correct for these noise sources in blood oxygenation level dependent (BOLD)-weighted time series. A dual‐echo sequence was used to simultaneously acquire data at both a short (TE=3.3ms) and a BOLD-weighted (TE=35ms) echo time. This approach is effectively “free,” using dead-time in the pulse sequence to collect an additional echo without affecting overall scan time or temporal resolution. The proposed correction method uses voxelwise regression of the short TE data from the BOLD-weighted data to remove noise variance. In addition to a typical resting state scan, non-compliant behavior associated with patient groups was simulated via increased head motion or physiological fluctuations in 10 subjects. Short TE data showed significant correlations with the traditional motion-related and physiological noise regressors used in current connectivity analyses. Following traditional preprocessing, the extent of significant additional variance explained by the short TE data regressors was significantly correlated with the average head motion across the scan in the resting data (r2=0.93, p<0.0001). The reduction in data variance following the inclusion of short TE regressors was also correlated with scan head motion (r2=0.48, p=0.027). Task-related data were used to demonstrate the effects of the short TE correction on BOLD activation time series with known temporal structure; the size and strength of the activation were significantly decreased, but it is not clear whether this reflects BOLD contamination in the short TE data or correlated changes in blood volume. Finally, functional connectivity maps of the default mode network were constructed using a seed correlation approach. The effects of short TE correction and low-pass filtering on the resulting correlations maps were compared. Results suggest that short TE correction more accurately differentiates artifactual correlations from the correlations of interest in conditions of amplified noise

Potential pitfalls when denoising resting state fMRI data using nuisance regression

In resting state fMRI, it is necessary to remove signal variance associated with noise sources, leaving cleaned fMRI time-series that more accurately reflect the underlying intrinsic brain fluctuations of interest. This is commonly achieved through nuisance regression, in which the fit is calculated of a noise model of head motion and physiological processes to the fMRI data in a General Linear Model, and the “cleaned” residuals of this fit are used in further analysis. We examine the statistical assumptions and requirements of the General Linear Model, and whether these are met during nuisance regression of resting state fMRI data. Using toy examples and real data we show how pre-whitening, temporal filtering and temporal shifting of regressors impact model fit. Based on our own observations, existing literature, and statistical theory, we make the following recommendations when employing nuisance regression: pre-whitening should be applied to achieve valid statistical inference of the noise model fit parameters; temporal filtering should be incorporated into the noise model to best account for changes in degrees of freedom; temporal shifting of regressors, although merited, should be achieved via optimisation and validation of a single temporal shift. We encourage all readers to make simple, practical changes to their fMRI denoising pipeline, and to regularly assess the appropriateness of the noise model used. By negotiating the potential pitfalls described in this paper, and by clearly reporting the details of nuisance regression in future manuscripts, we hope that the field will achieve more accurate and precise noise models for cleaning the resting state fMRI time-series

Is fMRI “noise” really noise? Resting state nuisance regressors remove variance with network structure

Noise correction is a critical step towards accurate mapping of resting state BOLD fMRI connectivity. Noise sources related to head motion or physiology are typically modelled by nuisance regressors, and a generalised linear model is applied to regress out the associated signal variance. In this study, we use independent component analysis (ICA) to characterise the data variance typically discarded in this pre-processing stage in a cohort of 12 healthy volunteers. The signal variance removed by 24, 12, 6, or only 3 head motion parameters demonstrated network structure typically associated with functional connectivity, and certain networks were discernable in the variance extracted by as few as 2 physiologic regressors. Simulated nuisance regressors, unrelated to the true data noise, also removed variance with network structure, indicating that any group of regressors that randomly sample variance may remove highly structured “signal” as well as “noise.” Furthermore, to support this we demonstrate that random sampling of the original data variance continues to exhibit robust network structure, even when as few as 10% of the original volumes are considered. Finally, we examine the diminishing returns of increasing the number of nuisance regressors used in pre-processing, showing that excessive use of motion regressors may do little better than chance in removing variance within a functional network. It remains an open challenge to understand the balance between the benefits and confounds of noise correction using nuisance regressors

Vascular physiology drives functional brain networks

We present the first evidence for vascular regulation driving fMRI signals in specific functional brain networks. Using concurrent neuronal and vascular stimuli, we collected 30 BOLD fMRI datasets in 10 healthy individuals: a working memory task, flashing checkerboard stimulus, and CO2 inhalation challenge were delivered in concurrent but orthogonal paradigms. The resulting imaging data were averaged together and decomposed using independent component analysis, and three “neuronal networks” were identified as demonstrating maximum temporal correlation with the neuronal stimulus paradigms: Default Mode Network, Task Positive Network, and Visual Network. For each of these, we observed a second network component with high spatial overlap. Using dual regression in the original 30 datasets, we extracted the time-series associated with these network pairs and calculated the percent of variance explained by the neuronal or vascular stimuli using a normalized R2 parameter. In each pairing, one network was dominated by the appropriate neuronal stimulus, and the other was dominated by the vascular stimulus as represented by the end-tidal CO2 time-series recorded in each scan. We acquired a second dataset in 8 of the original participants, where no CO2 challenge was delivered and CO2 levels fluctuated naturally with breathing variations. Although splitting of functional networks was not robust in these data, performing dual regression with the network maps from the original analysis in this new dataset successfully replicated our observations. Thus, in addition to responding to localized metabolic changes, the brain’s vasculature may be regulated in a coordinated manner that mimics (and potentially supports) specific functional brain networks. Multi-modal imaging and advances in fMRI acquisition and analysis could facilitate further study of the dual nature of functional brain networks. It will be critical to understand network-specific vascular function, and the behavior of a coupled vascular-neural network, in future studies of brain pathology

Arterial CO2 fluctuations modulate neuronal rhythmicity: Implications for MEG and fMRI studies of resting-state networks

A fast emerging technique for studying human resting state networks (RSNs) is based on spontaneous temporal fluctuations in neuronal oscillatory power, as measured by magnetoencephalography. However, it has been demonstrated recently that this power is sensitive to modulations in arterial CO2 concentration. Arterial CO2 can be modulated by natural fluctuations in breathing pattern, as might typically occur during the acquisition of an RSN experiment. Here, we demonstrate for the first time the fine-scale dependence of neuronal oscillatory power on arterial CO2 concentration, showing that reductions in alpha, beta, and gamma power are observed with even very mild levels of hypercapnia (increased arterial CO2). We use a graded hypercapnia paradigm and participant feedback to rule out a sensory cause, suggesting a predominantly physiological origin. Furthermore, we demonstrate that natural fluctuations in arterial CO2, without administration of inspired CO2, are of a sufficient level to influence neuronal oscillatory power significantly in the delta-, alpha-, beta-, and gamma-frequency bands. A more thorough understanding of the relationship between physiological factors and cortical rhythmicity is required. In light of these findings, existing results, paradigms, and analysis techniques for the study of resting-state brain data should be revisited

Multiparametric measurement of cerebral physiology using calibrated fMRI

The ultimate goal of calibrated fMRI is the quantitative imaging of oxygen metabolism (CMRO2), and this has been the focus of numerous methods and approaches. However, one underappreciated aspect of this quest is that in the drive to measure CMRO2, many other physiological parameters of interest are often acquired along the way. This can significantly increase the value of the dataset, providing greater information that is clinically relevant, or detail that can disambiguate the cause of signal variations. This can also be somewhat of a double-edged sword: calibrated fMRI experiments combine multiple parameters into a physiological model that requires multiple steps, thereby providing more opportunity for error propagation and increasing the noise and error of the final derived values. As with all measurements, there is a trade-off between imaging time, spatial resolution, coverage, and accuracy. In this review, we provide a brief overview of the benefits and pitfalls of extracting multiparametric measurements of cerebral physiology through calibrated fMRI experiments

Changes in arterial cerebral blood volume during lower body negative pressure measured with MRI

Cerebral Autoregulation (CA), defined as the ability of the cerebral vasculature to maintain stable levels of blood flow despite changes in systemic blood pressure, is a critical factor in neurophysiological health. Magnetic resonance imaging (MRI) is a powerful technique for investigating cerebrovascular function, offering high spatial resolution and wide fields of view (FOV), yet it is relatively underutilized as a tool for assessment of CA. The aim of this study was to demonstrate the potential of using MRI to measure changes in cerebrovascular resistance in response to lower body negative pressure (LBNP). A Pulsed Arterial Spin Labeling (PASL) approach with short inversion times (TI) was used to estimate cerebral arterial blood volume (CBVa) in eight healthy subjects at baseline and -40mmHg LBNP. We estimated group mean CBVa values of 3.13±1.00 and 2.70±0.38 for baseline and lbnp respectively, which were the result of a differential change in CBVa during -40 mmHg LBNP that was dependent on baseline CBVa. These data suggest that the PASL CBVa estimates are sensitive to the complex cerebrovascular response that occurs during the moderate orthostatic challenge delivered by LBNP, which we speculatively propose may involve differential changes in vascular tone within different segments of the arterial vasculature. These novel data provide invaluable insight into the mechanisms that regulate perfusion of the brain, and establishes the use of MRI as a tool for studying CA in more detail

Lag-Optimized Blood Oxygenation Level Dependent Cerebrovascular Reactivity Estimates Derived From Breathing Task Data Have a Stronger Relationship With Baseline Cerebral Blood Flow

Published: 15 June 2022Cerebrovascular reactivity (CVR), an important indicator of cerebrovascular health, is commonly studied with the Blood Oxygenation Level Dependent functional MRI (BOLD-fMRI) response to a vasoactive stimulus. Theoretical and empirical evidence suggests that baseline cerebral blood flow (CBF) modulates BOLD signal amplitude and may influence BOLD-CVR estimates. We address how acquisition and modeling choices affect the relationship between baseline cerebral blood flow (bCBF) and BOLD-CVR: whether BOLD-CVR is modeled with the inclusion of a breathing task, and whether BOLD-CVR amplitudes are optimized for hemodynamic lag effects. We assessed between-subject correlations of average GM values and within-subject spatial correlations across cortical regions. Our results suggest that a breathing task addition to a resting-state acquisition, alongside lag-optimization within BOLD-CVR modeling, can improve BOLD-CVR correlations with bCBF, both between- and within-subjects, likely because these CVR estimates are more physiologically accurate. We report positive correlations between bCBF and BOLD-CVR, both between- and within-subjects. The physiological explanation of this positive correlation is unclear; research with larger samples and tightly controlled vasoactive stimuli is needed. Insights into what drives variability in BOLD-CVR measurements and related measurements of cerebrovascular function are particularly relevant when interpreting results in populations with altered vascular and/or metabolic baselines or impaired cerebrovascular reserve.This work was supported by the Center for Translational Imaging at Northwestern University. The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research was supported by the Eunice Kennedy Shriver National Institute of Child Health and Human Development of the National Institutes of Health [K12HD073945]. KZ was supported by an NIH-funded training program [T32EB025766]. SM was supported by the European Union’s Horizon 2020 research and innovation program [Marie Skłodowska-Curie grant agreement No. 713673] and a fellowship from La Caixa Foundation [ID 100010434, fellowship code LCF/BQ/IN17/11620063]. CC-G was supported by the Spanish Ministry of Economy and Competitiveness [Ramon y Cajal Fellowship, RYC2017-21845], the Basque Government [BERC 2018-2021 and PIBA_2019_104], and the Spanish Ministry of Science, Innovation and Universities [MICINN; PID2019- 105520GB-100]

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function

Available online 24 June 2021.Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO 2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemo- dynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural ac- tivity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting- state data segments, and in data segments which added a 2–3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO 2 pressure: a breath-hold task to induce hypercapnia (CO 2 increase) and a cued deep breathing task to induce hypocapnia (CO 2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO 2 by systematically shifting the CO 2 regressor in time to optimize the model fit. This optimization inher- ently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO 2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemody- namic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.This research was supported by the Eunice Kennedy Shriver Na- tional Institute of Child Health and Human Development of the Na- tional Institutes of Health under award number K12HD073945. The pediatric dataset and cerebral palsy dataset were collected with sup- port of National Institutes of Health award R03 HD094615–01A1. The authors would like to acknowledge Marie Wasielewski and Carson Ingo for their support in acquiring these data. K.Z. was supported by an NIH-funded training program (T32EB025766). S.M. was supported by the European Union’s Horizon 2020 research and innovation pro- gram (Marie Sk ł odowska-Curie grant agreement No. 713673), a fel- lowship from La Caixa Foundation (ID 100010434, fellowship code LCF/BQ/IN17/11620063) and C.C.G was supported by the Spanish Ministry of Economy and Competitiveness (Ramon y Cajal Fellowship, RYC-2017- 21845), the Basque Government (BERC 2018–2021 and PIBA_2019_104) and the Spanish Ministry of Science, Innovation and Universities (MICINN; PID2019–105520GB-100)

The absolute CBF response to activation is preserved during elevated perfusion: Implications for neurovascular coupling measures

Functional magnetic resonance imaging (fMRI) techniques in which the blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) response to a neural stimulus are measured, can be used to estimate the fractional increase in the cerebral metabolic rate of oxygen consumption (CMRO2) that accompanies evoked neural activity. A measure of neurovascular coupling is obtained from the ratio of fractional CBF and CMRO2 responses, defined as n, with the implicit assumption that relative rather than absolute changes in CBF and CMRO2 adequately characterise the flow-metabolism response to neural activity. The coupling parameter n is important in terms of its effect on the BOLD response, and as potential insight into the flow-metabolism relationship in both normal and pathological brain function. In 10 healthy human subjects, BOLD and CBF responses were measured to test the effect of baseline perfusion (modulated by a hypercapnia challenge) on the coupling parameter n during graded visual stimulation. A dual-echo pulsed arterial spin labelling (PASL) sequence provided absolute quantification of CBF in baseline and active states as well as relative BOLD signal changes, which were used to estimate CMRO2 responses to the graded visual stimulus. The absolute CBF response to the visual stimuli were constant across different baseline CBF levels, meaning the fractional CBF responses were reduced at the hyperperfused baseline state. For the graded visual stimuli, values of n were significantly reduced during hypercapnia induced hyperperfusion. Assuming the evoked neural responses to the visual stimuli are the same for both baseline CBF states, this result has implications for fMRI studies that aim to measure neurovascular coupling using relative changes in CBF. The coupling parameter n is sensitive to baseline CBF, which would confound its interpretation in fMRI studies where there may be significant differences in baseline perfusion between groups. The absolute change in CBF, as opposed to the change relative to baseline, may more closely match the underlying increase in neural activity in response to a stimulus